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Abstract. We review the present status of global analyses of neutrino oscillations, 
taking into account the most recent neutrino data including the latest KamLAND 
and K2K updates presented at Neutrino2004, as well as state-of-the-art solar and 
atmospheric neutrino flux calculations. We give the two-neutrino solar + KamLAND 
results, and the two-neutrino atmospheric + K2K oscillation regions, discussing in 
each case the robustness of the oscillation interpretation against departures from the 
Standard Solar Model and the possible existence of non-standard neutrino physics. 
Furthermore, we give the best fit values and allowed ranges of the three-flavour 
oscillation parameters from the current worlds' global neutrino data sample and discuss 
in detail the status of the small parameters a = Am 2 OL /Am 2 TM as well as sin 2 #13, 
which characterize the strength of CP violating effects in neutrino oscillations. We also 
update the degree of rejection of four-neutrino interpretations of the LSND evidence 
in view of the most recent developments. 
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1. Introduction 



The discovery of neutrino masses by the combination of a variety of data from 



solar 0, H, i, 
and accelerator 




o masses by tne combination oi a 
1 1 Of ] , atmospheric HE 
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reactor 



13] neutrino experiments was the major recent achievement in 
astroparticle, high energy, and nuclear physics, which culminates a heroic effort dating 
back to about four decades. This has now firmly established the incompleteness of 
the Standard Model of electroweak interactions, expected on theoretical grounds since 
long ago 
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The determination of neutrino oscillation parameters 
is now a flourishing industry 23j which has finally entered the high precision age, 
with many experiments underway and a new generation coming. Apart from a careful 
understanding of solar and atmospheric neutrino fluxes, nuclear physics, neutrino cross 
sections and experimental response functions, the interpretation of the data relies heavily 
on the proper description of neutrino propagation properties both in the Sun and the 
Earth, including the so-called matter effects 24], |25| . 

The avalanche of data in a field where these have been traditionally so scarce, 
given the feebleness of neutrino interactions, has prompted a rush of phenomenological 
papers on the interpretation of neutrino data. A number of reviews are already in 
the market 
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Here we revisit the latest global analysis of neutrino 
oscillation parameters presented in Ref. 30[] in view of the most recent solar 211] and 
atmospheric 32J neutrino flux calculations and of the new data presented at Neutrino 
2004. 

This paper is organized as follows. In Sec. [2] we briefly discuss the analysis of the 
atmospheric neutrino data including the state-of-the-art three-dimensional calculation 
of the atmospheric neutrino flux given in Ref. 32| . In order to obtain the allowed ranges 
for the oscillation parameters Am^ TM and # ATM we combine the atmospheric neutrino data 
from the Super-K experiment with the accelerator data from the K2K experiment, which 
provides the first independent confirmation of the oscillation evidence from atmospheric 
neutrinos. We give also a short description of robustness of the atmospheric data with 
respect to non-standard neutrino interactions in Sec. 12.41 

In Sec. |3] we briefly describe the solar neutrino experiments and the KamLAND 
reactor neutrino experiment, and discuss the two- neutrino interpretation of the data 
described by Am| 0L and 9 SOL . We comment on the implications of the recent update of 



the Standard Solar Model (SSM) neutrino fluxes from Ref. (211]. Furthermore, in Sec. 13.41 
we discuss the robustness of the oscillation interpretation of solar neutrino data when 
going beyond the SSM framework or invoking non-standard neutrino properties. As 
examples we consider the effects of solar radiative-zone density fluctuations, convective- 
zone magnetic fields and the prospects of probing electromagnetic neutrino properties 
with current and future experiments. 

After discussing the dominant oscillations in the two-flavour approximation in 
Sees. [2] and |3] we devote Sec. |4] to the global three-neutrino analysis of the data, 
combining all current neutrino oscillation data except LSND. In addition to presenting 
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the best fit values and allowed ranges for the oscillation parameters, we update the 
status of the three-flavour parameters a = Am^/Am^, and #13 in view of the new 
data reported at Neutrino 2004. In Sec. [5] we give an update of the status of attempts to 
account for the LSND data in terms of four-neutrino oscillations. Finally we conclude 
and summarize our results in Sec. EJ A detailed discussion of the recent KamLAND 



results is given in Appendix A 



2. Leading oscillations with Am 



atm 



2.1. Atmospheric neutrino oscillations 

In 1998 the Super-K experiment obtained evidence for neutrino oscillations from the 
observation of the zenith angle dependence of their /i-like atmospheric neutrino data. 
This effect has been confirmed also by other atmospheric neutrino experiments, see e.g. 



Refs. [l2|, |l3j. Recently, Super-K has reported a dip in the L/E distribution of the 
atmospheric survival probability j^, which provides a clear signature for neutrino 
oscillations. 

In our atmospheric neutrino analysis we include the most recent charged- current 
atmospheric neutrino data from Super-K j34j, including the e-like and /i-like data 
samples of sub- and multi-GeV contained events (each grouped into 10 bins in zenith 
angle), as well as the stopping (5 angular bins) and through-going (10 angular bins) up- 
going muon data events. As previously, we do not use the information on v T appearance, 
multi-ring \i and neutral-current events since an efficient Monte-Carlo simulation of 
these data would require a more detailed knowledge of the Super-K experiment, and 
in particular of the way the neutral-current signal is extracted from the data. For 
details of our analysis see Refs. 35|, |36] and references therein. With respect to our 
previous atmospheric neutrino analysis 30j we have now taken into account the new 



three-dimensional atmospheric neutrino fluxes given in Ref. 32]. Furthermore, we 



have updated our statistical analysis following closely Ref. [371 ] . taking special care of 
systematical errors, like uncertainties in the neutrino fluxes and detection cross sections. 
Details on the \ 2 analysis can also be found in Ref. 38] . With these updates our results 
are in excellent agreement with the ones of the Super-K collaboration [34j . 

In Fig. [1] we show the results of our analysis of atmospheric data in the framework 
of two-flavour — ■> v T oscillations. The regions delimited by the hollow contours 
correspond to the allowed regions for the oscillation parameters sin 2 9 A 
The current best fit point occurs at 



and Am 2 TM . 



sin 2 9 ATM = 0.50 , Am 2 TM = 2.0 x 10~ 3 eV 2 (ATM data 



2 

ATM 



The main difference to the results of Ref. 30] are the relatively lower values of Am 
implied by the use of the new three-dimensional atmospheric neutrino fluxes reported 
in Ref. 32] instead of the one-dimensional Bartol fluxes 39] used previously. 
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Figure 1. Allowed (sin 2 9 A 



.,) regions at 90%, 95%, 99%, and 3cr C.L. for 2 



d.o.f. The regions delimited by the lines correspond to atmospheric data only, while 
for the colored regions also K2K data are added. The best fit point of atmospheric 
(atmospheric + K2K) data is marked by a triangle (star). Also shown is the A% 2 as a 
function of sin 2 9 ATM and Ato 2 tm , minimized with respect to the undisplayed parameter. 



2.2. The K2K accelerator experiment 



The KEK to Kamioka (K2K) long-baseline neutrino oscillation experiment 16j probes 
the disappearance oscillation channel in the same region of Am 2 as explored by 
atmospheric neutrinos. The neutrino beam is produced by a 12 GeV proton beam from 
the KEK proton synchrotron, and consists of 98% muon neutrinos with a mean energy of 
1.3 GeV. The beam is controlled by a near detector 300 m away from the proton target. 
Information on neutrino oscillations is obtained by the comparison of this near detector 
data with the content of the beam observed by the Super-Kamiokande detector at a 
distance of 250 km. 

The data sample called K2K-I [16| has been collected in the period from June 
1999 to July 2001 (4.8 x 10 19 protons on target). K2K-I has observed 56 events in 
Super-K, whereas 80. ll^ have been expected in the case of no oscillations. This 
gives a clear evidence for disappearance: the probability that the observed flux 
at Super-K is explained by a statistical fluctuation without neutrino oscillations is less 



than 1% 16j. Recently, at the Neutrino2004 conference new data from the K2K-II 
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period have been presented [17] . K2K-II started in fall 2002, and the released data 
corresponds to 4.1 x 10 19 protons on target, comparable to the K2K-I sample. From 
the combined analysis of K2K-I and K2K-II 108 events have been observed in Super-K, 
whereas 150.9^ have been expected for no oscillations. Out of the 108 events 56 are 
so-called single-ring muon events. This data sample contains mainly muon events from 
the quasi-elastic scattering v^+p — > fi+n, and the reconstructed energy is closely related 
to the true neutrino energy. Hence these data can be used for a spectral analysis. Using 
the Kolmogorov-Smirnov test the K2K collaboration finds that the observed spectrum is 
consistent with the spectrum expected for no oscillation only at a probability of 0.11%, 
whereas the spectrum predicted by the best fit oscillation parameters has a probability 
of 52% 0. 

In our re-analysis of K2K data we use the energy spectrum of the 56 single-ring 
muon events from K2K-I + K2K-II0. Similar as in Ref. [38J] we use a phenomenological 



parameterization for the spectrum expected for no oscillation. Adopting reasonable 
assumptions on the energy resolution function and on systematical errors, we fit the 
data divided into 15 bins in reconstructed neutrino energy, as given in Ref. [13]. The 
allowed regions for the oscillation parameters from K2K data are shown in Fig. |2] in 
comparison to the ones from atmospheric neutrino data. This figure illustrates that the 
neutrino mass-squared difference indicated by the disappearance observed in K2K is 
in perfect agreement with atmospheric neutrino oscillations. Hence, K2K data provide 
the first confirmation of oscillations with Am^ TM from a man-made neutrino source. 
K2K gives a rather weak constraint on the mixing angle due to low statistics in the 
current data sample. 

2.3. Atmospheric and K2K combined 

In Fig. [1] the allowed regions in the (sin 2 9 ATM , Am 2 TM ) plane from the combined analysis 
of K2K and Super-K atmospheric neutrino data are shown as shaded regions. As 
expected from Fig. [2] we find that, apart from providing an independent confirmation, 
K2K data start already to constrain the allowed region of Am 2 TM , whereas the 
determination of the mixing angle is completely dominated by atmospheric data. From 
the projections of the \ 2 onto the Am 2 TM and sin 2 9 ATM axes shown in Fig. [1] we find the 
best fit point 

sin 2 6 ATM = 0.5 , Am 2 TM = 2.2 x 10~ 3 eV 2 (ATM+K2K data) (2) 

with the corresponding allowed ranges at 3a (5a) for 1 d.o.f.: 

0.34 (0.27) < sin 2 ATM < 0.66 (0.73) , (3) 
1.4 (0.85) x 10~ 3 eV 2 < Am 2 TM < 3.3 (4.2) x 10~ 3 eV 2 . 

Note that despite the downward shift of the atmospheric mass-splitting implied by 



the new neutrino fluxes from Ref. [32] our new quoted value for Amf TM in Eq. ([2]) is 



| We cannot use the full K2K data sample of 108 events, since not enough information is available to 
analyze these data outside the K2K collaboration. 
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Figure 2. Allowed K2K regions in the (sin 9 atm, Aml TU ) plane at 90%, 95%, 99%, and 
3cr C.L. for 2 d.o.f. The hollow lines delimit the region determined from the atmospheric 
data only. The star (triangle) corresponds to the K2K (atmospheric) best fit point. 



statistically compatible both with our previous result |30j and the value obtained by 
the new Super-K L/E analysis 



33 



Let us remark that the K2K constraint on Am^ TH 
from below is important for future long-baseline experiments, since the performance of 
such experiments is drastically affected if AmL, were in the lower part of the 3a range 



allowed by atmospheric data (see e.g. Ref. |40j). 



2.4- Robustness of oscillation interpretation: atmospheric neutrinos 

We now turn to the issue of the robustness of the oscillation interpretation of the 
atmospheric neutrino data. Non-standard physics may in principle affect atmospheric 
neutrino fluxes, as well as neutrino propagation and detection cross sections 26 



Apart from the issue of their theoretical viability [41| neutrino decays have been 
considered since long ago j42|. Non-standard interactions arising from the SU(2] <S> U(l) 
charged and neutral currents 



20 



or from new particles [43j, |44], |45|, |46|, |47|, |48|, |49|, as 
well as quantum mechanical decoherence 38] might also affect atmospheric neutrino 
results. Althou gh strongly rejected by recent atmospheric data as the dominant 



mechanisms 33j, |38|, |50j, non-standard phenomena might still be present at a sub- 
leading level in addition to oscillations and, to this extent, have some impact on the 
determination of the oscillation parameters. 

In the following we illustrate the stability of the measurement of AmJ TM and 
sin 2 # ATM by assuming the presence of non-standard interactions of the neutrinos with 
earth matter. New neutrino interactions beyond the Standard Model are a natural 



Status of global fits to neutrino oscillations 




Figure 3. Effective flavour-changing NSI operator. 



feature in most neutrino mass models 51] and can be of two types: flavour-changing 
(FC) and non- universal (NU). These interactions (called NSI for short) may be 
schematically represented as effective dimension-6 terms of the type eGf, as illustrated 
in Fig. 131 where e specifies their sub-weak strength. Such interactions may arise from a 
nontrivial structure of charged and neutral current weak interactions characterized by a 
non- unitary lepton mixing matrix [20| . These gauge-induced NSI may lead to flavour and 

Alternatively, 



45 



47 



48 



49 



CP violation, even with massless or degenerate neutrinos 

such non-standard neutrino interactions may also arise in models where neutrino masses 
are "calculable" from radiative corrections 43, 



44 



and in some supersymmetric models 
with broken R parity j52|. Finally, in supersymmetric unified models, the strength of 
non-standard neutrino interactions may be a calculable renormalization effect [461 ] . 

The impact of non-standard neutrino interactions on atmospheric neutrinos was 
considered in Ref. 50| treating the NSI strengths as free phenomenological parameters. 
This analysis takes into account both the effect of — > v T oscillations (OSC) as well 
as the existence of non-standard neutrino-matter interactions (NSI) in this channel. In 
addition to the standard term in the Hamiltonian describing oscillations an term -£/nsi is 
introduced, accounting for an effective potential induced by the NSI with earth matter: 



H 



NSI 



±V2G F N f 



e 



e e' 



(4) 



Here +(— ) holds for neutrinos (anti-neutrinos) and e and e' parameterize the NSI: 
\f2GpNfE is the forward scattering amplitude for the FC process + / — > v T + / and 
V2GpNfe' represents the difference between v^+f and v T +f elastic forward scattering. 
The quantity Nf is the number density of the fermion / along the neutrino path. For 
definiteness we take for / the down-type quark. 

In Fig. H] we show the results of a fit to atmospheric neutrino data for the four 
parameters Am 2 TM , sin 2 20 AXM , e, and e'%. In the left panels the pure oscillation case 
(e = e' = 0) is compared to the case where also some NSI are allowed. We find that 
the x 2 improves slightly (2.4 units) in the presence of NSI, but the determination of 
the oscillation parameters is practically unaffected. This is an important result, since 
it shows that the allowed ranges derived for AmJ TM and sin 2 29 ATM are rather stable 
with respect to non-standard physics [jf In turn, the high preference of the data for 



Here we assume that the parameter e is real; for the more general case of complex s, see Ref. [53 
Note that the analysis of Ref. [HB] is based on the neutrino fluxes of Ref. , and in addition to the 
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Figure 4. Atmospheric neutrino oscillation parameters from [50j fitted in the 
presence of non-standard interactions using the Bartol flux • We show the behaviour 
of the x 2 as a function of the four parameters Am^ TM (upper left), sin 2 20 ATM (lower 
left), the FC parameter e (upper right), and the NU parameter e' (lower right). In 
each panel the x 2 is minimized with respect to the three undisplayed parameters. In 
the left panels we show also the pure oscillation case (OSC), with e and e' set to zero. 



oscillations allows to set strong bounds on NSI. From the right panels of Fig. H] we 
deduce the bounds at 3cr 



-0.03 <e< 0.02, |er'| < 0.05 . 



(5) 



Before closing we stress that, apart from its intrinsic theoretical importance, the 
study of non-standard neutrino interactions has an astrophysical interest, as these can 
affect the propagation of neutrinos in a variety of astrophysical environments, such 



54 1551 ] and pulsars [56(]. They could lead to "deep- inside" neutrino 



as supernovae 

conversions, in addition to those expected from conventional neutrino oscillations. 



Super-K data described in Sec. 12.11 also the up-going muon data from the MACRO experiment [12| is 
used. This explains the slight difference between the best fit values of Am 2 TM from Fig. 0] for the pure 
oscillation case and Eq. (Tl]). 
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3. Leading oscillations with Am^ 

3.1. Solar neutrino oscillations 

In this section we consider oscillations of solar neutrinos in the two-flavour framework. 
A detailed discussion of experimental and theoretical aspects of solar neutrino physics 



can be found in other contributions to this volume [57|, l58j. Therefore, we focus in 
the following on the determination of the oscillation parameters from the fit to solar 
neutrino data. In our analysis we take into account the rates of the chlorine experiment 
at the Homestake mine [Jy] (2.56±0.16±0.16 SNU), the most up-to-date results 59] of 
the gallium experiments SAGE 0,0 (66.9 ±f% SNU) and GALLEX/GNO jjH,0 
(69.3 ± 4.1 ± 3.6 SNU), as well as the 1496-day Super-K data sample 0] in the form of 
44 bins (8 energy bins, 6 of which are further divided into 7 zenith angle bins). From the 
SNO experiment we include the most recent data from the salt phase [Io| in the form 
of the neutral current (NC), charged current (CC) and elastic scattering (ES) fluxes, as 
well as the 2002 spectral day/night data 0, 0] (17 energy bins for each day and night 
period). 



The analysis methods used here are similar to the ones described in Refs. 30|, |36J] and 
references therein, including the use of the so-called pull approach for the \ 2 calculation, 



as described in Ref. 6j|. In this method all systematic uncertainties are included by 
introducing new parameters in the fit and adding a penalty function to the x 2 . For 
example, for each of the eight solar neutrino fluxes a parameter is introduced, and the 
predictions from the Standard Solar Model including the correlated errors are taken into 



account by means of a penalty function. The method described in Ref. [6l| is extended 

as 



in two respects. First, it is generalized to the case of correlated statistical errors [62 
necessary to treat the SNO-salt data. Second, we do not consider the x 2 on ly U P to first 
order in the pulls, but instead each pull parameter is treated exactly to all orders. This 
is particularly interesting in the case of the solar 8 B flux. In our approach it is possible 
to include the SSM prediction for this flux as well as the SNO NC measurement on the 
same footing, without pre-selecting a particular value, as implied by expanding around 
the predicted value. In this way the fit itself can choose the best compromise between 
the SNO NC data and the SSM prediction. 

In Fig. [5] we compare the allowed regions for the oscillation parameters using the 



solar neutrino fluxes given in the BP00 SSM 63| and the recent update (BP04) presented 



in Ref. |31j. One finds that the change in the flux predictions has a negligible impact 
on the allowed regions. This illustrates that thanks to the good experimental accuracy 
the determination of the oscillation parameters is rather robust with respect to changes 
in the SSM. The current best fit values for solar neutrino oscillation parameters are 

sin 2 6 SOL = 0.29 , Am 2 0L = 6.0 x 10~ 5 eV 2 (solar data, BP04).(6) 



Also the rejection against maximal solar mixing is 5.6a, the same as found previously [30] 
using the BP00 solar model. This is the significance at which bi-maximal models of 
neutrino mass, such as the CP conserving version of the neutrino unification model 
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Figure 5. Allowed regions from all solar neutrino data at 90%, 95%, 99%, and 3cr 
C.L. for 2 d.o.f. in the plane of sin 2 6 W and Am 2 0L . The regions delimited by the 
curves correspond to the BP00 SSM 63|, whereas for the colored regions the BP04 
SSM 31[ has been used. Also shown is Ax 2 as a function of sin 2 9 S0L and Am 2 OL , 
minimized with respect to the undisplayed parameter. The labeled contours denote 
constant CC/NC ratio in the SNO experiment. 



given in Ref. 



64 



are ruled out. Despite the fact that the BP04 model does not imply 
any significant change in the solar neutrino parameters, we will use it in our subsequent 
discussion of solar neutrino results. 



3.2. The KamLAND reactor neutrino experiment 

The KamLAND experiment is a reactor neutrino experiment with its detector 

located at the Kamiokande site. Most of the v e flux incident at KamLAND comes 
from nuclear plants at distances of 80 — 350 km from the detector, making the average 
baseline of about 180 kilometers, long enough to provide a sensitive probe of the LMA 
solution of the solar neutrino problem. The KamLAND collaboration has for the first 
time measured the disappearance of neutrinos traveling to a detector from a power 
reactor. They observe a strong evidence for the disappearance of neutrinos during their 
flight over such distances, giving the first terrestrial confirmation of the solar neutrino 
anomaly and also establishing the oscillation hypothesis with man-produced neutrinos. 
In KamLAND the reactor anti-neutrinos are observed by the process u e + p — > 
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e + + n, where the delayed coincidence of the prompt energy from the positron and 
a characteristic gamma from the neutron capture allows an efficient reduction of 
backgrounds. The neutrino energy is related to the prompt energy by E v = E pT +A—m e , 
where A is the neutron-proton mass difference and m e is the positron mass. In the lower 
part of the energy spectrum there is a relevant contribution from geo-neutrino events 



to the signal (see, e.g., Refs. |65|, |66j). To avoid large uncertainties associated with the 



geo-neutrino flux an energy cut at 2.6 MeV prompt energy is applied for the oscillation 
analysis. 

First results from KamLAND were published in Ref. 14|. In the period from March 
to October 2002 data corresponding to a 162 ton-year exposure have been collected, and 
after all cuts 54 anti-neutrino events remained in the final sample. This number has to 
be compared with 86.8 ± 5.6 reactor neutrino events predicted for no oscillations and 
0.95 ± 0.99 background events, which gives a probability that the KamLAND result 



is consistent with the no-disappearance hypothesis of less than 0.05% [14] . Recently 



new results have been presented by KamLAND 15]. With a somewhat larger fiducial 



volume of the detector an exposure corresponding to 766.3 ton-year has been obtained 
between March 2002 and January 2004 (including a reanalysis of the 2002 data from 
Ref. 14j). A total of 258 events has been observed, in comparison to the expectation 
of 356.2 ± 23.7 reactor neutrino events in the case of no disappearance and 7.5 ± 1.3 
background events. This leads to a confidence level of 99.995% for v e disappearance, 
and the averaged survival probability is 0.686 ± 0.044(stat) ± 0.045(syst). Moreover 
evidence for spectral distortion consistent with oscillations is obtained [151 ]. 



It was shown in Ref. 



67] in relation with the first KamLAND data that most 
information on the neutrino oscillation parameters sin 2 # S0L and Am 2 OL can be extracted 
from the data by using an event-based likelihood analysis. This requires the knowledge 
of the prompt energy of each observed event. Unlike to the first data sample it is not 
possible to recover this information for the latest KamLAND data sample from publicly 
available material. Therefore only x 2 -analyses based on binned data can be performed 
outside the KamLAND collaboration. Details of our KamLAND simulation based on 
2002 data can be found in Refs. |67|, [68|, a discussion of our updated analysis method 
for the current data sample is given in Appendix A Here we briefly summarize the 



main features: Instead of the traditional bins of equal size in E pr we use the data 



binned equally in 1/E pi , which can be obtained from Ref. [15], and which allows to 



extract more relevant information on the oscillation parameters. We adopt the improved 



anti- neutrino flux parameterization from Ref. 69], and include various systematic 
errors associated to the neutrino fluxes, reactor fuel composition and individual reactor 
powers. Furthermore, we include matter effects, a careful treatment of backgrounds 
and information on the average contribution to the total reactor neutrino signal as a 
function of the distance to the detector. 

The KamLAND allowed regions for sin 2 ^sol and Am 2 0L are shown in Fig. [6] in 
comparison to the regions from solar data. One observes beautiful agreement between 
KamLAND data and the region implied by the LMA solution to the solar neutrino 
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Figure 6. Allowed regions from KamLAND data at 90%, 95%, 99%, and 3ct C.L. for 
2 d.o.f. for the 2002 data sample [13] (left panel) and the latest 766.3 ton-year data 
sample [li} (right panel). The regions delimited by the lines correspond to solar data. 
The KamLAND best fit points are marked with a star, the solar best fit point with a 
dot. 



problem, which in this way has been singled out as the only viable one. Before this 
experiment we had a very complex pattern of alternative oscillation solutions like LOW, 



SMA, or VAC, see e.g., Refs. |35|, [36]. All of these are completely ruled out by the 
KamLAND data. From this point of view the KamLAND experiment has played a 
fundamental role in the resolution of the solar neutrino problem. 

Comparing the left and right panel of Fig. [6] one can appreciate the improvement 
implied by the recent KamLAND data. The allowed region is drastically reduced, and 
in particular the mass-squared difference is very well determined by recent data. We 
find the best fit point 

Am s 2 OL = 8.1 +0 i x 1(T 5 eV 2 (la) _ . , A , _ 

" OL ~ - 3 v 1 (KamLAND data). (7) 

sin 2 fl SOL = 0.29 1 ; 1 ; 

The solution at low values of Am 2 0L is present at the 99% C.L., with the local minimum 
at Am 2 OL = 1.6 x 1(T 5 eV 2 and sin 2 SOL = 0.31 and A^ 2 = 7.5 relative to the best fit 
point. The so-called high-LMA solution is present only marginally at the 3cr level at 
Am 2 OL = 1.7 x lCr 4 eV 2 , sin 2 9 SOh = 0.25 and A^ 2 = 11.3. Although not statistically 
significant, we note that the small matter effect favours slightly values of sin 2 9 SOL < 0.5 
over the "mirror solution" with sin 2 9 SOL > 0.5. The high values of Am 2 OL > 3 x 10~ 4 eV 2 



in the averaging regime, which have been allowed by first KamLAND data 14j at 



90% C.L. are now ruled out at more than 3a. This indicates that the spectral distortion 
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Figure 7. Allowed regions from combined solar and KamLAND data at 90%, 95%, 
99%, and 3cr C.L. for 2 d.o.f.. The regions shown with lines correspond to solar data 
only, from Fig. 



associated to neutrino oscillations has been observed [15(] (see also Appendix A for 
further discussion). 



3.3. Solar and KamLAND combined 

Under the fundamental assumption of CPT invariance we can directly compare the 
information obtained from solar neutrino experiments and the KamLAND reactor 
experiment. In Fig. [7] we show the allowed regions from the combined solar and 
KamLAND data. The current best fit point of the global analysis occurs at 

sin 2 # SOL = 0.29 , Am 2 OL = 8.1x 1(T 5 eV 2 (Solar+KamLAND data). (8) 

From the projections of the \ 2 onto the Am^ OL and sin 2 # S0L axes also shown in Fig. [7] 
we find the following allowed ranges at 3a (5<r) for 1 d.o.f.: 

0.23 (0.19) < sin 2 6 SOL < 0.37 (0.45) , (9) 
7.3 (6.7) x 10~ 5 eV 2 < Am 2 OL < 9.1 (9.9) x 10~ 5 eV 2 . 

As expected from Fig. [6] the determination of the mixing angle is completely dominated 
by the solar neutrino data, whereas KamLAND significantly reduces the allowed range 
for Am 2 OL . Also the so-called high-LMA region, which previously was present around 




Figure 8. Effect of random matter density fluctuations with a correlation length of 
Lq = 100 km on the electron-neutrino survival probability for LMA oscillations. 



Am 2 SOL ~ 1.4 x 10~ 4 eV 2 at 3cx (see e.g., Ref. [30]) is now ruled out with a Ax 2 = 22.1 
with respect to the global minimum, which corresponds to an exclusion at about 4.3a 

An 



30, 



36 



for 2 d.o.f.. Notice also that the day/night data are treated as previously 
improved analysis of the day /night asymmetry data along the lines followed by the 
Super-K collaboration in Ref. 70| would lead to an even more pronounced rejection 
of the high Am^ region. Currently not enough information is available to reproduce 
this result outside the Super-K collaboration. Finally we note that the use of the BP04 
fluxes has not changed significantly the allowed regions. 



3-4- Robustness of oscillation interpretation: solar neutrinos 

The oscillation interpretation of solar neutrino data depends both on astrophysical input 
(the model of the Sun) as well as on the physics characterizing both the propagation 
as well as neutrino interaction properties. Either may differ from the Standard Solar 
Model and Standard Electroweak Model expectations. How robust is the oscillation 
interpretation of solar neutrino data in view of this? 



3.4-1- Beyond the Standard Solar Model In the following we will briefly discuss 
consequences of departures from the Standard Solar Model. The effect of varying 
solar neutrino fluxes has been widely discussed, thus, study we consider the 

possibility of solar density fluctuations. Even though this possibility was suggested in 



a number of papers 7l|, |72j, |73J] it has been traditionally neglected for several reasons. 



First, helioseismic measurements constrain deviations of solar properties from Standard 
Solar Model predictions at better than the percent level. Second, preliminary studies 



of the implications for neutrino oscillations of radiative-zone helioseismic waves 73 



indicated that they were unlikely to have observable effects. Third, no other known 
sources of fluctuations seemed to have the properties required to influence neutrino 
oscillations. 

Recently all of these points have been re-examined, with the result that the 
presence of solar fluctuations appears more likely than previously thought. First, direct 
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Figure 9. Solar neutrino oscillation parameters with an arbitrary noise amplitude and 
a correlation length Lq = 100 km from combining solar data with 2002 KamLAND 
data (left panel) and 2004 KamLAND data (right panel). The contour lines refer to 
solar data only. 



helioseismic bounds turn out to be insensitive to fluctuations whose size is as small as 
those to which neutrinos are sensitive 0, [tHI] typically several hundreds of km. Second, 
recent studies have shown how such solar density fluctuations can arise near the solar 
equatorial plane in the presence of magnetic fields deep within the solar radiative zone 



due to a resonance between Alfven waves and helioseismic g-modes 76 



It has been shown in Ref. 23] that such density fluctuations can affect neutrino 
propagation in an important way. The effect of random matter density fluctuations on 
the electron-neutrino survival probability for LMA oscillations has been shown to be 
sizable if the correlation length L (we take L Q = 100 km) is comparable to the neutrino 
oscillation length in the Sun. This is illustrated in Fig. [SJ The fluctuation's amplitude 
£ at the position of neutrino resonance is zero in the left panel, and is £ = 4% and 
£ = 8% in the middle and right panels, respectively. The corresponding solar neutrino 
oscillation parameters obtained in our global fit are shown in Fig. [9] before and after 
the new KamLAND data presented at Neutrino 2004. One sees that these new data 
have a rather strong impact on the stability of the oscillation parameters. With 2002 
KamLAND data an additional allowed region was present for Am s 2 OL ~ 2 x 10- 5 eV 2 at 
the 99% C.L. 78] . whereas the new KamLAND data pin down the oscillation parameters 
such that the allowed regions in Fig. [9] are practically stable, and the new solar noise- 
induced solution appears only marginally at the 3cr confidence level. 

Conversely, as shown in Ref. [78], the quality of current solar neutrino measurements 
after the SNO-salt and KamLAND results is sufficiently good as to place important 
constraints on fluctuations in the solar medium deep within the solar radiative zone. 
In other words, neutrinos may be used as an astrophysical probe of the solar interior, 
beyond the framework of the Standard Solar Model. As illustrated in Fig. [10] density 




Figure 10. Bounds on random matter density fluctuations for a correlation length of 
Lo = 100 km from solar neutrino data, solar + KamLAND data, and solar data for 
oscillation parameters fixed at the best fit point. 



fluctuations are strongly constrained if the correlation length lies in the range of several 
hundred km. Comparing the curves for free and fixed oscillation parameters one notes 
that the bounds on fluctuations have already become rather stable, a situation that 
may still improve when the accuracy on Am;? OL and 6 S0L gets improved by future 
neutrino experiments. Because oscillations are sensitive to correlation lengths which 
are so short, such solar neutrino results will complement the constraints that come from 
helioseismology. 



3.4-2. Beyond solar neutrino oscillations: Spin Flavour Precession In extensions of 
the Standard Model (SM) neutrino masses are in general accompanied also by non- 
standard interactions and/or electromagnetic properties. In the minimal extension of 
the Standard Model with Dirac neutrino masses one expects very tiny neutrino magnetic 



moments (MMs) 79], well below current experimental sensitivities. However, the 
theoretically preferred case of Majorana neutrinos leads to potentially larger transition 
magnetic moments closer to the present sensitivities. These can affect neutrino 
propagation properties in the Sun beyond the oscillation mechanism, due to the possible 
presence of solar magnetic fields. Alternatively, they can affect the determination of 
neutrino oscillation parameters due to non-standard neutrino cross sections inside the 
detectors. 

The most general form of the electromagnetic current of massive (Majorana) 
neutrinos has been given in Ref. 



80 



The magnetic piece is characterized by a 3 x 3 
complex antisymmetric matrix, the so-called Majorana transition moment (TM) matrix, 
that contains MMs as well as electric dipole moments of the neutrinos. Their existence 
would affect neutrino propagation inside the solar convective zone due to a spin-flavour 
precession (SFP) effect 80|, |8lJ, |82[. This in general depends on the assumed magnetic 
field profile. In order to quantify the extent with which the oscillation regions can be 
altered by the sub- leading spin-flavour precession effect, a \ 2 analysis was performed 
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in Ref. [83] taking into account the global solar + KamLAND disappearance data. We 
assumed that neutrino conversions are driven mainly by LMA oscillations, and used 
the same self-consistent [84j convective-zone solar magnetic field profile employed in 



Ref. 85] . The results we obtain indicate that, even though small, current bounds 
on neutrino magnetic moments and solar magnetic fields still leave room for slight 
modifications in the determinations of solar neutrino oscillation parameters, in the 
presence of large magnetic moments 0. 

However, in the general Majorana case, where theory may give rise to higher 
moments, there is a characteristic feature of the spin flavour precession which will 
lead to more stringent constraints, and hence increase the robustness of the oscillation 
parameter determination. The argument is based on the presence of anti-neutrinos in 
the solar flux 80|, |8l|, |82j. Recently the KamLAND collaboration .86] has reported a 
result which greatly improves the bound on an anti-neutrino component in the solar 
flux from 0.1% of the solar boron v e flux to 2.8 x 1(T 2 % at the 90% C.L., about 30 



times better than the recent Super-K limit 87]. This implies that, in practical terms, 
extremely good stability of the solar neutrino oscillation parameters against the possible 
existence of sub- leading SFP conversions is obtained. As a result, we conclude that 
solar neutrinos oscillation parameters can be inferred without any reference to neutrino 



83 for 



magnetic properties nor solar magnetic fields. We refer the reader to Ref. 
quantitative details. 

All in all, our analysis of solar neutrino data implies indicates pretty good stability 
of the oscillation parameter determination for the case of Majorana neutrinos, due to 
the solar anti-neutrino limit from KamLAND. In contrast for the special case of Dirac 
neutrinos this limit does not apply and the determination of oscillation parameters 
is potentially more fragile. However, we note that the gauge theoretic expectations for 
Dirac magnetic moments are typically lower than those for Majorana neutrino transition 
moments. 



3.4-3. Beyond SM neutrino cross sections: Constraining neutrino magnetic moments 
Neutrino transition magnetic moments are basic properties of neutrinos 80|. Although 
they do not substantially affect neutrino propagation, even in the presence of solar 
magnetic fields, non-trivial electromagnetic neutrino properties could still show up in the 
detection process and to this extent affect the determination of oscillation parameters. 
Experiments based on the neutrino detection via neutrino-electron elastic scattering 
are a sensitive probe of the electromagnetic properties. In Ref. 88] it was shown that 
current data from solar neutrinos (in particular from Super-K) in combination with 
reactor neutrino-electron scattering data provides strong bounds on all the elements of 
the TM matrix (for similar analyses see Ref. 89_|, |9CJ, |9jJ, [92] ) . 



In several experiments such as Super-K, Borexino and some reactor experiments [93 



94], 125], neutrinos are detected via the elastic neutrino-electron scattering, whose 



% Note that our analysis of solar neutrino data applies also to the special case of Dirac neutrinos. 
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electromagnetic cross section is [96|, 197 
da Prn a 2 n ( 1 1 



alT mln% \T E v ' ^ ' (10) 

where // e g is an effective MM [98] , T denotes the kinetic energy of the recoil electron and 
E v is the incoming neutrino energy. The electromagnetic cross section adds to the weak 
cross section and allows to extract information on the TM matrix, which we denote by 
A in the following. Taking into account the antisymmetry of A for Majorana neutrinos, 
it is useful to define vectors A by Xj k = £ju^-h where Xjk are the elements of the TM 
matrix in the mass basis. The effective MM square // 2 ff takes on different forms in the 
cases of solar and reactor neutrino experiments. The detailed derivation of the following 
expressions can be found in Ref. 88]. For the case of solar neutrino experiments one 
obtains the effective MM square 

^ma = I A| 2 - |A 2 | 2 + P% (|A 2 | 2 - |Ai| 2 ) , (11) 

where P 2 ^ corresponds to the probability that an electron neutrino produced in the core 
of the sun arrives at the detector as the mass eigenstate V\ in a two-neutrino scheme. 
In contrast, the yU 2 ff relevant in reactor experiments is given as 

/4 = |A| 2 -cos 2 fl SOL |Ai| 2 - sin 2 6> SOL |A 2 | 2 - sin20 SOL |Ai||A 2 | costf, (12) 

where 5 = arg(A^A 2 ) is the relative phase between A x and A 2 . 

In the following we discuss the constraints on neutrino TMs from solar and reactor 
neutrino experiments [88) . The \ 2 obtained from the data is minimized with respect to 
all TM parameters except the modulus |A|. To take into account the physical boundary 
|A| > we use Bayesian methods to calculate an upper bound on |A|. Let us stress 
that these bounds apply to all elements of the TM matrix, including MMs and electric 
dipole moments of all neutrino flavours, since |A| 2 = |Ai| 2 + 1 7V2 1 2 + | VV3 1 2 . Furthermore, 
since |A| is independent of the basis these bounds apply also for the TMs in the flavour 
basis. 

In the left panel of Fig. [11] we show contours of the 90% C.L. bound on |A| in the 
(tan 2 9 SOL , Am SOL ) plane for the combination of solar and reactor data. We note that in 
the upper parts of the LMA region, the solar data alone give already a strong bound on 
|A|, see Ref. 88] for details. In contrast, for low Am 2 0L values the inclusion of reactor 
data plays an important role in improving the bound. From our analysis we find at 90% 
C.L.: 



. 1 3.4 x IO^Vb (solar + KamLAND data) 
1 ' 1.7x10 _1 Vb (solar + KamLAND + reactor data), 



where for each value of |A| we have minimized the x 2 with respect to tan 2 6 l SOL and 
Am 2 



"SOL- 



Finally we note that the Borexino experiment [99(] will improve these bounds 
roughly by one order of magnitude. This experiment is mainly sensitive to the solar 7 Be 
neutrino flux, which will be measured by elastic neutrino-electron scattering. Therefore, 
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tan 2 e tan 2 G 

Figure 11. Contours of the 90% C.L. bound on |A| in units of 10 _10 /i-B from combined 
solar and reactor data (left panel) and after 3 years of Borexino data-taking (right 
panel). The gray (light) shaded region is the 3cr LMA region obtained in the global 
analysis of solar neutrino data (best fit point marked with a triangle), whereas the green 
(dark) one corresponds to the 3er region obtained after including the KamLAND results 
(best fit point marked with a star). The dashed line in the right panel corresponds to 
P e i = 0.5 for 7 Be neutrinos, and shows the strongest attainable limit. 



Borexino is similar to Super-K, the main difference is the mono-energetic line of the 
7 Be neutrinos, with an energy of 0.862 MeV, which is roughly one order of magnitude 
smaller than the energies of the 8 B neutrino flux relevant in Super-K. Thanks to the 
lower neutrino energy the sensitivity to electromagnetic properties is increased, as can 



be seen from Eq. (110p . Details about our Borexino simulation can be found in Ref. 88 
At the best fit point one finds the sensitivity 

|A| < 0.29 x IO^Vb at 90% C.L. (14) 

after three years of Borexino data taking. In the right panel of Fig. [11] we show contours 
of the 90% C.L. bound in the (tan 2 8 SOh , Amg 0L ) plane. 



4. Three— flavour neutrino oscillations 



4-1. Global three-neutrino analysis 



In this section the three-neutrino oscillation parameters are determined from a global 
analysis of the mos t rece nt neutrino oscillation data. For earlier three-neutrino analyses 
see Refs. |35|, llOCH . 1 1 11 ] . To fix the notation, we define the neutrino mass-squared 
differences AmL, = Am?, = mi, — m? and AmL„ = Am?n = mi — m\, and use the 



SOL 



convenient form of the parameterization for the leptonic mixing matrix given in Ref. [20 
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parameter 


best fit 


2a 


3a 


4a 


Am§! [10" 5 eV 2 ] 


8.1 


7.5-8.7 


7.2-9.1 


7.0-9.4 


Amli [10- 3 eV 2 ] 


2.2 


1.7-2.9 


1.4-3.3 


1.1-3.7 


sin 2 


0.30 


0.25-0.34 


0.23-0.38 


0.21-0.41 


sin 2 23 


0.50 


0.38-0.64 


0.34-0.68 


0.30-0.72 


sin 2 0i3 


0.000 


< 0.028 


< 0.047 


< 0.068 



Table 1. Best-fit values, 2cr, 3cr, and 4cr intervals (1 d.o.f.) for the three-flavour 
neutrino oscillation parameters from global data including solar, atmospheric, reactor 
(KamLAND and CHOOZ) and accelerator (K2K) experiments. 




and now adopted as standard by the PDG 102 

' C13C12 S12C13 

U = — S12C23 — S23S13C12 C23C12 — S23S13S12 S23C13 J . (15) 
^ S23S12 - S13C23C12 -S23C12 - S13S12C23 

where Cy = cos0jj and Sij = sin0jj. Furthermore, we use the notations 812 = S ol an d 
023 = 0atm- Because of the hierarchy Am 2 0L <C Am 2 TM it is agood approximation to 
set Am 2 OL = in the analysis of atmospheric and K2K dataQ and to set Am 2 TM to 
infinity for the analysis of solar and KamLAND data. This implies furthermore that 
the effect of apossible Dirac CP-violating phase 20j in the lepton mixing matrix can 
be neglected Q We perform a general fit to the global data in the five- dimensional 
parameter space s 2 2 , sl 3 , sf 3 , Am|i, Am| 1; and show projections onto various one- or 
two-dimensional sub-spaces. 

We include in our analysis the global solar neutrino oscillation data from all solar 
neutrino experiments and the KamLAND reactor experiment as described in Sec. [31 
the atmospheric neutrino data from Super-K, as well as spectral data from the K2K 
long-baseline experiment (see Sec. [2]). In additio n we take into account in our fit the 



105 



constraints from the CHOOZ reactor experiment 

The results of the global three-neutrino analysis are summarized in Fig. [12] and 
in Tab. [TJ In the upper panels of the figure the A% 2 is shown as a function of 
the parameters sin 2 0i2, sin 2 023, sin 2 0i3, Am 2 1; Am^, minimized with respect to the 
undisplayed parameters. The lower panels show two-dimensional projections of the 
allowed regions in the five-dimensional parameter space. The best fit values and 
the allowed ranges of the oscillation parameters from the global data are given in 
Tab. [TJ This table summarizes the current status of the three-flavour neutrino oscillation 
parameters. 



+ See Ref. [1 031 ] for a two-mass scale analysis of atmospheric data. 

* The two Majorana phases [20J do not show up in oscillations but do appear in lepton number violating 
processes [104 1. 
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Figure 12. Projections of the allowed regions from the global oscillation data at 90%, 
95%, 99%, and 3<r C.L. for 2 d.o.f. for various parameter combinations. Also shown is 
A% 2 as a function of the oscillation parameters sin 2 #12, sin 2 $23, sin 2 #13, Am^, Am^, 
minimized with respect to all undisplayed parameters. 



4-2. The small parameters a 



Aml 0L /Am 2 ATM and 6 13 



Genuine three-flavour effects are associated to the mass hierarchy parameter a = 
^ m soL/^ m ATM an d the mixing angle #13 . In particular, in a three-neutrino sc heme 



20 



106 



CP violation disappears in the limit where two neutrinos become degenerate 
and in the limit where #13 — > 0. We discuss in this subsection the present status of these 
small parameters. 

In Fig. [13]the Ax 2 from the global data is shown as a function of the mass hierarchy 
parameter a. Also shown in this figure is the Ax 2 as a function of the parameter 
combination asin2# 12 , since to leading order in the long baselin e v e 



oscillation 



probability solar parameters appear in this particular combination 1071 . Il08| . We obtain 
the following best fit values and 3a intervals: 



a = 0.035 
asin26 , 12 = 



0.032 



0.024 < a < 0.060, 

0.022 < asin2# 12 < 0.054. 



(16) 



Let us now discuss the status of the mixing angle #13, which at the moment is 
the last unknown angle in the three-neutrino leptonic mixi ng m atrix. Only an u pper 
bound exists, which used to be dominated by the CHOOZ 1051 ] and Palo Verde [109 



reactor experiments. Currently a large effort is put to determine this angle in future 





Figure 14. Ax 2 profiles projected onto the sin 2 6*13 axis, for solar, KamLAND, 
solar+KamLAND, atmospheric+K2K+CHOOZ, and for the global data. 



experiments (see, e.g., Refs. [261. 1271. l40j). 



In Fig. [T3]we show the A% 2 as a function of sin 2 #13 for different data sample choices. 
From this figure we find the following bounds at 90% C.L. (3<r) for 1 d.o.f.: 

{ 0.041 (0.079) (solar+KamLAND) 
0.029 (0.067) (CHOOZ+atmospheric+K2K) (17) 
0.022 (0.047) (global data) 

We find that the new data from KamLAND have a surprisingly strong impact on 
this bound. Before the 2004 KamLAND data the bound on sin 2 # 13 from global data 
was dominated by the CHOOZ reactor experiment, together with the determination of 
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Figure 15. Upper bound on sin 2 13 (1 d.o.f.) from solar+KamLAND+CHOOZ 
data as a function of Am 2 ml . The dashed (solid) curve corresponds to the 90% (3cr) 
C.L. bound, the thin curves have been obtained with 2002 KamLAND data, whereas 
the thick curves follow from the recent 2004 KamLAND update. The light (dark) 
shaded region is excluded by CHOOZ data alone at 90% (3er) C.L. The horizontal 
line corresponds to the best fit value of Am 2 TM from atmospheric + K2K data given 
in Eq. ([2|, and the hatched regions are excluded by atmospheric + K2K data at 3a 
according to Eq. (J3J. 



Am^ from atmospheric data (see e.g., Ref. 30|). However, using most recent data the 
combined bound from solar+KamLAND becomes comparable to the CHOOZ bound, 
and these data contribute notably to the final bound. A detailed discussion of the 
reason for such improvement on sin 2 #13 from the 2004 KamLAND data is given in 
Appendix A.2| One reason is the rather strong signal for spectral distortion in the 



current sample. 

As noted in Ref. 30j the bound from solar+KamLAND is especially important for 
the relatively lower values of AmJ TM implied by the use of the new three-dimensional 
atmospheric fluxes 32j, since the CHOOZ bound on sin 2 #13 deteriorates quickly when 
Am 2 TM decreases (see Fig. H5j) . Such loosening in sensitivity for low Am 2 TM values is 
prevented first, by the lower bound on Am 2 TM from K2K (see Fig. [2]) and second, by 
the bound from solar+KamLAND, which is independent of Am 2 TM . In Fig. [151 we show 
the upper bound on sin $13 as a function of Am 2 TM from CHOOZ data alone compared 
to the bound from an analysis including solar and reactor neutrino data (CHOOZ and 
KamLAND). One finds that, although for larger Am 2 TM values the bound on sin 2 # 13 is 
dominated by the CHOOZ data, for Am 2 TM < 2 x 10~ 3 eV 2 the solar + KamLAND data 
start being important. For illustration we show in Fig. [15] also the bound as implied 
by the old 2002 KamLAND data to highlight the improvement of the new data. We 
note that, as before, the bound from solar data is rather stable under the recent SSM 
update. 
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Figure 16. Sections of the three-dimensional allowed regions in the (sin 2 9 SOL , Am 2 0L ) 
plane at 90%, 95%, 99% and 3c C.L. for 3 d.o.f. for various sin 2 6\z values from solar 
data (lines) and solar+KamLAND data (colored regions). The local minima in each 
plane from solar+KamLAND (solar only) data are marked by filled (open) dots. 



Let us discuss in some more detail the constraint on sin 2 (9 13 from solar data, which 
emerges from a subtle interplay of various solar neutrino observables. In Fig.[16]we show 
the results of a three parameter fit (sin 2 9 S0L , Am 2 0L , sin 2 #13) to solar and KamLAND 
data. Allowed regions are shown for various values of sin 2 6*13 in the (sin 2 9 S0L , Am 2 0L ) 
plane with respect to the global minimum. Note that here we calculate the allowed 
regions at a given confidence level for 3 d.o.f. The shape of A% 2 from solar data shown 
in Fig. [H] can be understood from Fig. [T6j Indeed one observes that for solar data 
increasing #13 can be compensated to some extent by increasing Am 2 OL . Since solar 
data disfavours large values of Am 2 0L the bound improves. Also the combination with 
KamLAND has a similar effect, since recent KamLAND data essentially fix Am 2 OL at 
roughly 8 x eV 2 such that the continuous rise of Am 2 OL with sin 2 #13 preferred by 
solar data is prevented. On the other hand solar data breaks a correlation of sin 2 #13 and 
sin 2 9 12 in the KamLAND data (see Appendix A. 2 ), which again leads to an improvement 
of the combined bound on sin 2 # 13 . 

The difference in the day / night solar neutrino fluxes due to the regenerati on effect in 



110 



111 



the earth in the three-flavour framework has been considered recently in Refs. 
This observable may provide valuable infor matio n on #13 in the context of future solar 



neutrino experiments like UNO or Hyper-K 112 



5. Four— neutrino oscillations and LSND 

In addition to the strong evidence for oscillations due to the mass-squared differences 
Am 2 0L and Am 2 TM there is also a hint for o scillations with a much larger mass-squared 



difference from the LSND experiment [113j. This accelerator experiment performed at 
Los Alamos observed 87.9 ±22.4 ±6.0 excess events in the — > v e appearance channel, 
corresponding to a transition probability of P = (0.264 ± 0.067 ± 0.045)%, which is 
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Figure 17. The two classes of six four-neutrino mass spectra, (3+1) and (2+2). 



3.3a away from zero. To explain this signal with neutrino oscillations requires a 



mass-squared difference Ami 



1 eV . Such a value is inconsistent with the mass- 



squared differences required by solar/KamLAND and atmospheric/K2K experiments 
within the standard three-flavour framework. In this section we consider four-neutrino 
schemes, where a sterile neutrino 
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1161 is added to the three active ones to 



provide the additional mass scale needed to reconcile the LSND evidence. We include 
in our anal ysis data f rom the LSN D experiment, as well as from short-baseline (SBL) 
accelerator 1171 . Ill8l ] an d rea ctor 1051 . 1 1091 . Ill9j ] experiments reporting no evidence 
for oscillations (see Ref. 12CH ] for details of ou r SBL da ta analysis). We update our 
previous four-neutrino analyses (see, e.g., Refs. [12lL Il22| ) by including the most recent 
solar and KamLAND 15] data, the improved atmospheric neutrino fluxes 32] and latest 
data from the K2K long-baseline experiment (l7| . 



5.1. A common parameterization for four-neutrino schemes 

Four-neutrino mass schemes are usually divided into the two classes (3+1) and (2+2), as 
illustrated in Fig. [T7J We note that (3+1) mass spectra include the three-active neutrino 
scenario as limiting case. In this case solar and atmospheric neutrino oscillations are 
explained by active neutrino oscillations, with mass-squared differences A?rig 0L and 
Am1 TM , and the fourth neutrino state gets completely decoupled. We will refer to such 
limiting scenario as (3+0). In contrast, the (2+2) spectrum is intrinsically different, as 
the sterile neutrino must take part in either solar or in atmospheric neutrino oscillations, 
or in both. 

Neglecting CP violation, neutrino oscillations in four-neutrino schemes are 
generally described by 9 parameters: 3 mass-squared differences and 6 mixing a ngle s in 



123 



m 



the lepton mixing matrix [20j. We use the parameterization introduced in Ref. 
terms of Am;? 0L , 9 S0L , AmJ TU , 9 ATM , AmJ SND , 6*lsnd- These 6 parameters are similar to 
the two-neutrino mass-squared differences and mixing angles and are directly related to 
the oscillations in solar, atmospheric and the LSND experiments. For the remaining 3 
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parameters we use r) s ,r) e and d^. These quantities are defined by 



d r 



\ U Q i\ 2 with i E solar mass states, 



with % G atmospheric mass states, 



(19) 



where a = e, /x, r, s. Note that in (2+2) schemes the relation r] a = d a holds, whereas 
in (3+1) r] a and d a are independent. The physical meaning of these parameters is the 
following: r\ a is the fraction of v a participating in solar oscillations, and (1 — d a ) is 
the f racti on of v a participating in oscillations with Am 2 ATM (for further discussions see 
Ref. 123]). For the analysis we adopt the following approximations: 

(i) We make use of the hierarchy Am;? OL <C Am^ TM <C AmJ SND . This means that for 
each data set we consider only one mass-squared difference, the other two are set 
either to zero or to infinity. 

(ii) In the analyses of solar and atmospheric data (but not for SBL data) we s et r\ e = 1, 



whi ch is justified because of strong constraints from reactor experiments 1051 . Il09 
119 . 



Within this approximation the parameter structure of the four-neutrino analysis 
gets rather simple. The parameter dependence of the four data sets solar, atmospheric, 
LSND and NEV is illustrated in Fig. [TSJ In this section, except where explicitly noted 
otherwise, we tacitly consider KamLAND as part of the solar data sample and K2K as 
part of the atm ospheric data sample. The NEV data set contains the experi ment s 
KARMEN [mi, CDHS |nl|, Bugey jllj, CHOOZ [l05t and Palo Verde [l09 



reporting no evidence for oscillations. We see that only r] s links solar and atmospheric 
data and d^ links atmospheric and NEV data, while LSND and NEV data are coupled 
by AmL SND and 9 LSND . With the definitions f|T8|) and f fl9|) and in our approximation the 
parameter str uctu re shown in Fig. HH] holds for both types of mass spectra, (3+1) as 



well as (2+2) [123 
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Figure 19. Left: A\ 2 as a function of from solar data before the SNO salt- 
phase results, from current solar data, and from solar+KamLAND data. Middle: 
AxLl) Ax1 tm+ k2k+sbl and Xgiobai as a function of ij s in (2+2) oscillation schemes. 
The dashed line corresponds to atmospheric and K2K data only (without SBL data). 
Right: Axatm+k2k as a function of d^. 



5.2. (2+2): ruled out by solar and atmospheric data 

The strong preference for oscillations into active neutrinos in solar and atmospheric 
oscillations |36j leads to a direct conflict in (2+2) oscillation schemes. We will now show 
that thanks to recent solar neutrino data (in particular from the SNO-salt phase (lQ| ) 
in combination with the KamLAND experiment [14| . and the latest Super-K data on 



atmospheric neutrinos HI] the tension in the data has become so strong that (2+2) 



oscillation schemes are essentially ruled out 0. 

In the left panel of Fig. [19] we show the A% 2 from solar neutrino data as a function 
of rj s , the parameter describing the fraction of the sterile neutrino participating in solar 
neutrino oscillations. It is clear from the figure that the improved determination of 
the neutral current event rate from the solar 8 B flux implied by the salt enhanced 



measurement in SNO 10] substantially tightened the constraint on a sterile contribution: 
the 99% C.L. bound improves from from rj s < 0.44 for pre-SNO-salt to rj s < 0.31 
(BP00). The boron flux predicted in the current BP04 SSM is slightly larger than the 
NC flux measured in SNO, which leaves more room for a sterile component in the solar 
neutrino flux. Indeed, using the BP04 SSM the bound deteriorates slightly to r] s < 0.33 
at the 99% C.L. This effect illustrates that in schemes beyond minimal three-flavour 
oscillations the data still shows some sensitivity to the theoretical SSM input. Although 
KamLAND on its own is insensitive to a sterile neutrino contamination, it contributes 
indirectly to the bound because of the better determination of Am;? 0L and 9 S0L [36|. The 
combined analysis leads to the 99% C.L. bound 

7] s < 0.25 (solar + KamLAND, BP04). (20) 



f For an earlier four-neutrino analysis of solar and atmospheric data see Ref. 



124 1 
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In contrast, in (2+2) schemes atmospheric data prefer values of r] s close to 1. From 
the combined analysis of Super-K atmospheric data, K2K and SBL neutrino data we 
obtain the bound r] s > 0.75 at 99% C.L., in clear disagreement with the bound from 
solar data. In the middle panel of Fig. [19] we show the Ax 2 for solar data and for 
atmospheric+K2K combined with SBL data as a function of rj s . Note that the main 
effect comes from atmospheric+K2K data; SBL experiments contribute only marginally, 
as visible from the dashed line. From this figure we also see that the "solar+KamLAND" 
and the "atm+K2K+SBL" allowed domains overlap only at Xpc = 17.2, i.e. at the 4.1a 
level. 

In the middle panel of Fig. [19] we also show the "global" y 2 function defined as 
follows: 

X 2 (Vs) = AXsol+kam(^) + AxLm+k2k+sbl(^) • (21) 



In Refs. [12 lL Il25| we have proposed a statistical method to evaluate the disagreement 
of different data sets in global analyses. The parameter goodness of fit (PG) makes 
use of the x 2 defined in Eq. ( 12~TT) . This criterion is very useful to evaluate the 



GOF of the combination of data sets, avoiding dilution by the large num ber of data 



points, as it happens for the usual GOF criterion (for details see Ref. [1251 ] ) . We 
find Xpg = Xmin = 26.1, which corresponds to 5.1cr. We conclude that (2+2) mass 
schemes are ruled out by the disagreement between the latest solar and atmospheric 
neutrino data. This is a very robust result, independent of whether LSND is confirmed 
or disproved 0. 

Let us note that we now obtain a slightly smaller Xpg than previously 122] , and 



the disagreement gets slightly weaker. The reason is that in the present analysis we 



have not included the atmospheric neutrino data from the MACRO experiment 12J. As 



discussed in Ref. [36|, these data have some sensitivity to rj s and enhance the rejection 
against a sterile component in atmospheric oscillations. Since relevant information to 
perform a consistent combined analysis of Super-K and MACRO using the new fluxes 
is presently not available, we prefer to use only Super-K data, which allows us to derive 
simpler and more robust resultsjjjl Let us mention that we are also neglecting neutral 
current and tau appearance data from Super-K, which would also increase the rejection 

| Sub-leading effects beyond the approximations adopted here should not affect this result significantly. 
Allowing for add itional parameters to vary at the percent level might change the ratio of some 
observables |l26| |. however, we expect that the absolute number of events relevant for the fit will 
not change substantially. 

§ The reason for this is the presence of some tension between theoretical predictions from three- 
dimensional fluxes and MACRO as well as Super-K thru-going muon data (see discussion in second 
reference of [13]). This tension did not appear with older one-dimensional fluxes, and makes it 
particularly important to properly take into account all the correlations between different data sets. For 
Super-K detailed information on the correlation between various sub-samples (sub-GeV, multi-GeV, 
stopping and thru- going muons) can be extracted from Ref. [37j], whereas the relevant information to 
include MACRO is not available. Since in this analysis we are more interested in the robustness rather 
than in the strength of our results, we restrict our analysis to Super-K data only. A consistent way of 
combining Super-K and MACRO data using three-dimensional fluxes is however under investigation. 
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Figure 20. Upper bound on sin 2 29 LSND from NEV, atmospheric and K2K neutrino 
data in (3+1) schemes. The bound is calculated for each Atol SND using the A\ 2 for 
1 d.o.f. The dotted line corresponds to the bound at 99% C.L. without K2K and using 
one-dimensional atmospheric flux es. A lso shown are the regions allowed at 99% C.L. 
(2 d.o.f.) from global LSND data [HI and decay-at-rest (DAR) LSND data (l28 |. 



against a sterile component, since the detailed informations which are needed to use 
these data are presently not available outside the Super-K collaboration. 



5.3. (3+1): strongly disfavoured by SBL data 



It is known for a long time 120l . 1 1291 . Il30l . Il3ll . 1 1321 . Il33l . Il34j | that (3+ 1) mass schemes 
are disfavoured by the comparison of SBL disappearance data [1181 . Ill9| with the LSND 
result. The reason is that in (3+1) schemes the relation sin 2 2# LSND = 4d e dp holds, 
and the parameters d e and dp (see Eq. (jl9]) ) are strongly constrained by u e and Vp 
disa ppear ance experiments, leading to a double suppression of the LSND amplitude. In 
Ref. 132| it was realized that the up-down asymmetry observed in atmospheric [i events 
leads to an additional constraint on dp. The Ax 2 {dp) from the fit to atmospheric+K2K 
data is shown in the right panel of Fig. [191 We find that the use of the new atmospheric 
fluxes {32| as well as K2K data [13] considerably strengthen the constraint on dp-, the 
new bound dp < 0.065 at 99 % C.L . decreases roughly a factor 2 with respect to the 
previous bound dp < 0.13 0, fl27| . as implied by fitting atm ospheric data with one- 
dimensional fluxes 39] and without K2K. Following Ref. 1271 ] we show in Fig. [20] the 
upper bound on the LSND oscillation amplitude sin 2 2# LSND from the combined analysis 
of NEV and atmospheric neutrino data. This figure illustrates that the improvement 
implied by the stronger bound on dp is mostly relevant for lower values of Am 2 SND . From 
this figure we see that the bound is incompatible with the signal observed in LSND at the 
95% C.L. Only marginal overlap regions exist between the bound and global LSND data 
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SOL ATM LSND NEV 


Xpg 


PG 


(3+1) 
(2+2) 


0.0 0.4 5.7 10.9 
5.3 20.8 0.6 7.3 


17.0 
33.9 


1.9 x 10~ 3 (3.1a) 
7.8 x 10~ 7 (4.9a) 



Table 2. Parameter GOF and the contributions of different data sets to Xpg m (3+1) 
and (2+2) neutrino mass schemes. 



if both are taken at 99% C.L. Using only the decay-at-rest LSND data sample 1281 ] the 
disagreement is even more severe. These results show that (3+1) schemes are strongly 
disfavoured by SBL disappearance data. 



5.4- Comparing (3+1), (2+2) and (3+0) hypotheses 



Using the methods developed in Ref. 123| we perform a global fit to the oscillation 
data in the four-neutrino framework. This approach allows to statistically compare 
the different hypotheses. Let us first evalu ate the GOF of (3+1) and (2+2) spectra 



using the PG method described in Ref. 



125 



We divide the global oscillation data into 



the four data sets SOL, ATM, LSND and NEV. Then we use the PG method to test 
the statistical c omp atibility of these data sets assuming a given neutrino mass scheme. 



Following Ref. 



121 



we consider 



X" = Axsol^sol, Am 2 SOL ,r] s ) + Ax 2 ATM (0 ATM 

LSND) Am LSNB , dp, Ve) + AXlsnd(#lsnd, A^lsnd) > 

where A Xx = Xx ~ Ode)**. (X = SOL, ATM, NEV, LSND). In Tab. we show the 
contributions of the 4 data sets to Xpg = Xmm f° r (3+1) and (2+2) oscillation schemes. 
As expected we observe that in (3+1) schemes the main contribution comes from SBL 
data due to the tension between LSND and NEV data in these schemes. For (2+2) 
oscillation schemes a large part of Xpg comes from solar and atmospheric data, due to 
the rejection against a sterile neutrino contribution of these two data sets, as discussed 
in Sec. 15.21 The contribution fr om N EV data in (2+2) comes mainly from the tension 
between LSND and KARMEN 1281 ] . which does not depend on the mass scheme. 



^sol, Am 2 OL , rj s ) 



(22) 



The parameter goodness of fit is obtained by evaluating \pg f° r 4 d.o.f. |l25 



This number of degrees of freedom corresponds to the 4 parameters r) s , d^, LSND , Aml SND 
describing the coupling of the different data sets (see Eq. ( 1221) and Fig. H8"]l . The best 
GOF is obtained in the (3+1) case. However, even in this best case the PG is o nly . 19% . 



121 



122 



This PG value is slightly reduced with respect to our previous result 0.56% 
because of the improved limit on d^ from atmospheric+K2K data (see Sec. 15.31) . The 
PG of 7.8 x 10 -7 for (2+2) schemes shows that these mass schemes are essentially ruled 
out by the disagreement between the individual data sets. A s me ntioned in Sec. 15.21 the 
exclusion of (2+2) schemes is slightly weaker as previously 122] since we now do not 
include MACRO data in our analysis. 

Although we have seen that none of the four-neutrino mass schemes provides a 
good fit to the global oscillation data including LSND, it is interesting to consider 
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the relative status of the three hypotheses (3+1), (2+2) and the three-active neutrino 
scenario (3+0). This can be done by comparing the x 2 value of the best fit point - which 
occurs for the (3+1) scheme - with the ones corresponding to (2+2) and (3+0). First 
we observe that (2+2) schemes are strongly disfavoured with respect to (3+1) with a 
Ax 2 = 16.9. For 4 d.o.f. this is equivalent to an exclusion at 3.1a. Furthermore, we find 
that (3+0) is disfavoured with a A% 2 = 17.5 (corresponding to 3.2a for 4 d.o.f.) with 
respect to (3+1). This reflects the high statistical significance of the LSND result, since 
in a (3+0) scheme no effect is predicted for LSND. 

To summarize, we find that four-neutrino schemes do not provide a satisfactory fit 
to the global data. The strong rejection of non-active oscillation in the solar+KamLAND 
and atmospheric+K2K neutrino data rules out (2+2) schemes, irrespective of whether 
LSND is confirmed or not. Using an improved goodness of fit method especially sensitive 
to the combination of data sets we find that (2+2) schemes are ruled out at the 4.9a 
level. On the other hand (3+1) spectra are disfavoured by the disagreement of LSND 
with short-baseline disappearance data, leading to a marginal GOF of 1.9 x 10~ 3 (3.1a). 
Should LSND be confirmed it would be very desirable to have more data on v e and/or 
SBL disappearance to decide about the status of (3+1) schemes. In that case a positive 
signal is predicted right at the sensitivity edge of existing experiments. 

More drastic attempts to reconcile the LSND signal with the rest of neutrino 



oscillation data have been reviewed in Ref. [26|]. For example, in Ref . 135j | a five- 
neutrino scheme is invoked to reconcile all the data. In Ref. 136| it has been 
shown that even the rather drastic assumption of CPT violation in a three-neutrino 
framework [1371 1 1381 ] does not provide a satisfactory description of the global neutrino 
data set including LSND. Similarly, an interpretation o f the LSND signal in terms of a 



139 



non-standard muon decay is disfavoured by KARMEN 

We conclude that currently no convincing explanation for the LSND result exists, 
and it remains a puzzle how to reconcile this evidence with the rest of the data. It 
is therefore very important to settle this issu e exp erimentally. A confirmation of the 



LSND signal by the MiniBooNE experiment 140 
require some novel physics ideas. 



would be very exciting and would 



6. Summary and conclusions 

We have given a brief review of the status of global analyses of neutrino oscillations, 
taking into account the latest neutrino data, including the most recent updates of 
KamLAND and K2K presented at Neutrino2004, as well as state-of-the-art solar and 
atmospheric neutrino flux calculations. We presented two-neutrino solar + KamLAND 
results, as well as two-neutrino atmospheric + K2K oscillation regions, and a discussion 
in each case of the robustness with which the oscillation hypothesis can be established, 
in view of possible modifications. These might come from the assumed theoretical 
fluxes, the non- validity of the Standard Model neutrino interaction cross sections or the 
existence of non-trivial neutrino propagation properties beyond oscillations. As case 
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studies we have mentioned the robustness of the solar neutrino oscillation hypothesis 
vis a vis the possible existence of radiative-zone solar density fluctuations, nonzero 
convective-zone solar magnetic fields and neutrino transition magnetic moments. For 
the atmospheric + K2K analysis we have considered explicitly the robustness of the 
oscillation hypothesis against the possible existence of flavour or universality violating 
non-standard neutrino interactions. 

Furthermore, we have performed a fit to the most recent world neutrino data sample 
in the three-flavour framework. The results of this global analysis are summarized in 
Fig. [12] and Tab. [IJ where we give the best fit values and allowed ranges of the three- 
flavour oscillation parameters. In addition we discussed in detail the status of the small 
parameters a = Amg OL /Am| TM and sin 2 #13, which characterize the strength of CP 
violating effects in neutrino oscillations, highlighting the improvement of the bound on 
sin 2 #13 implied by the inclusion of the recent KamLAND data. Finally, we gave a review 
over the current status of four-neutrino interpretations of the LSND anomaly, in view 
of the most recent experimental and theoretical advances. 

All in all, we can say beyond reasonable doubt that neutrino masses, discovered 
through atmospheric neutrino oscillations, have now also been confirmed in the 
solar neutrino oscillation channel thanks to the important input of the KamLAND 
experiment. Theory-wise, while the SSM was necessary in order to establish the need 
for physics beyond the Standard Model, it has now been made to some extent irrelevant 
by the high precision of the experiments which currently dominate the determination of 
solar neutrino oscillation parameters. The next goal in the agenda is the determination 
of the small parameter sin 2 # 13 that characterizes the strength of CP violating effects in 
neutrino oscillations, and the exploration of the Majorana nature of the neutrino which 
will be sensitive to the other leptonic CP phases. 
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Appendix A. The new KamLAND results 



In this appendix we discuss in some detail the most recent data from KamLAND 15 



In |Appendix A. 1| we give the details of our data analysis, whereas in |Appendix A. 2 
discuss the physics results. 



wc 



Appendix A.l. Analysis details of the new data 



We have modified our previous KamLAND analysis methods 63, |68| applied to the first 
published KamLAND data 14] in various aspects to take into account the characteristics 
of the new data. First, to analyze KamLAND data one has to know the contribution 
of the various power reactors to the signal. We extract the relevant information from 
Fig. 1(b) of Ref. 15], where the no-oscillation signal is given as a function of the distance 
to the detector. Second, we use an improved parameterization of the anti-neutrino 
flux emitted by the isotopes 235 U, 239 Pu, 238 U and 241 Pu in the nuclear reactors. Third, 
we include 2.69 accidental background events in the lowest energy bin. For the 4.8 ±0.9 
background events expected from the beta-decay of 9 Li and 8 He we assume regular beta- 
spectra with endpoints of 13.6 and 10.7 MeV, respectively (see, e.g., Fig. 11 of Ref. [401]). 
Fourth, we include the small matter effects. 

In general most information can be extracted from data by un-binned Likelihood 
methods (see Ref. 67] for the case of the first KamLAND data). Unfortunately it is not 
possible to obtain event-based energy information for the current KamLAND sample, 
and one has to stick with binned data outside the collaboration. Traditionally data is 
given in bins of equal size in the prompt energy E pr (see Fig. 2(b) of Ref. 15]). However, 
it turns out that in the case of KamLAND more information can be obtained if data 
is binned equally in l/E pT . The relevant information can be extracted from Fig. 3 of 
Ref. 15] , where the ratio of the observed spectrum to the expectation for no oscillation 
is given in 13 bins of ^Okm/E^ [MeV]. We show the data binned in 1/E pr as well as 
in E pr in Fig. IA1I The 1/E pi binning is more useful for two reasons. First, it is more 
natural to make the bins smaller in the region of many events (low energy) and wider 
in the high energy region, where there are very few events. This maintains more energy 
information in the region of greater statistics. Second, since the frequency of neutrino 
oscillations is proportional to 1/E V this binning is more appropriate for the signal we 
are interested. 

Concerning the statistical analysis, we adopt a Poisson x 2 -function, and make 
extensive use of the pull-method to implement various systematical errors. In the 
overall normalization uncertainty we include only the detector-specific contributions 
by summing up the errors of the left column of Tab. 1 from Ref. 15], which gives 
c"det = 5.47%. The uncertainties associated to the anti- neutrino flux are treated 
according to the method presented in Ref. 69]. We include an uncertainty on the 
thermal power (2%) and the fuel composition (1%) of each individual reactor, as well as 



the spectral uncertainty of the emitted anti-neutrino fluxes. In agreement with Ref. 69 
we find that flux related uncertainties play only a minor role in the KamLAND analysis. 
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Figure Al. Data divided by the expectation for no oscillations for equal bins in 1/E pt 
(triangles) and in E pl (dots). Also shown is the spectrum for the best fit point (solid), 
the local best fit points in the low Am 2 (dash-dotted) and high Am 2 (dash-dot-dotted) 
regions, and for a large value of sin 2 613 (dashed). 



Appendix A. 2. Discussion of the KamLAND results 

In Fig. IA1I we show the predictions for the probabilities in the various bins Pi = 
^osc/^no-osc compared to the data P° hs = A^ bs /iV* _ 0SC . From this figure one can 
see that the probabilities are rather low for low energies, whereas in the region around 
4.5 MeV even slightly more events than expected for no oscillations have been observed. 
This is a very characteristic pattern indicating rather strong spectral distortion, and 
only the best fit parameters (solid line) can reproduce this shape. All other parameter 
choices shown in the figure have problems to accommodate the high data points in the 
middle of the spectrum. This holds for the solution around Am 2 ~ 1.6 x 10~ 5 eV 2 as 
well as for the high-LMA solution at Am 2 ~ 1.7 x 10~ 4 eV 2 (compare Fig. [6]). The 
rejection power to these two "solutions" is significantly increased by the 1/E pr binning 
with respect to the E pT binning: In our analysis we obtain a Ax 2 = 7.5(11.3) for the 
local minimum at 1.6(17) x 10~ 5 eV 2 , whereas with the traditional E pT binning we find 
Ax 2 = 4.0(8.5). In particular, to disfavour the high-LMA region the more precise energy 
information at low energies provided by the 1/E pi bins is crucial, whereas for the low 
region the high energy part seems to be important. 

Let us now discuss the bound on sin 2 # 13 from the KamLAND data. From Fig. IA1I 
one can see that turning on sin 2 #13 leads to a flatter energy spectrum and it gets more 
difficult to accommodate the low probabilities for low energies simultaneously with the 
high data points in the middle part of the spectrum. It is also clear from the relevant 




Figure A2. Left panel: x 2 from KamLAND data as a function of sin 2 #13 including 
the unphysical region sin 2 #13 < 0. Right panel: allowed regions at 90%, 95%, 99% 
and 3ct of sin 2 9 13 and sin 2 U (2 d.o.f.) for fixed Am 2 OL = 8.1 x 1CT 5 eV 2 from 
KamLAND data (dashed), solar data (solid), and KamLAND+solar data combined 
(shaded regions). 



three-flavour survival probability in vacuum 

Pee = 1 - ~ sin 2 26 13 - cos 4 13 sin 2 26 12 sin 2 ^f 1 ^ (A.l) 

that values of sin 2 # 13 > suppress the oscillatory term. In fact, if a fit to the KamLAND 
data is performed without imposing the constraint sin 2 # 13 > one finds a best fit point 
within the unphysical region at sin 2 # 13 = —0.13 (see Fig. IA2|) . Because of the rather 
large difference in the probability between the low and medium energy bins the fit is 
improved by a A% 2 ~ 3.5 by allowing an enhancement of the oscillatory term in Eq. ( 1A.1I) 
due to values of cos 4 6*i 3 > 1, i.e., sin 2 # 13 < 0. Fig. IA2I illustrates that this fact leads 
to a rather strong bound on sin 2 #i 3 if the analysis is restricted to the physical region. 
Just from the statistical power of the data one would expect a ler error on sin 2 #i 3 of 
a ~ 0.07. However, due to the particular fluctuation observed in the actual data the 
best fit point lies in the unphysical region. This implies that for sin 2 #i 3 = there is 
already some tension in the fit, and a significantly smaller error of a = 0.02 is obtained 
within the physical regionjjj] 

Finally, in the right panel of Fig. IA2I we illustrate, why the combination of 
KamLAND and solar data leads to a further notable improvement of the bound on 
sin 2 #i 3 . From Eq. ( 1A.1I) one expects for KamLAND a correlation between sin 2 #i 3 and 
sin 2 2# 12 . Solving Eq. (lA.lj) for sin 2 2# 12 and expanding up to first order in sin 2 # 13 one 



I We note that in such a case a reliable bound can be calculated by performing a Monte Carlo simulation 
of many synthetic data sets. This is beyond the scope of the present article, and we throughout use the 
naive method of calculating bounds by considering Ax 2 -values as implied by Gaussian statistics and 
restricting the analysis to the physical region. 



Status of global fits to neutrino oscillations 



36 



13 13 

without C background with C background 




0.2 0.4 0.6 0.8 10 0.2 0.4 0.6 0.8 1 
.2 .2 

sin 12 sin 12 

Figure Bl. Allowed (sin 2 6» SOL , Am£ OL ) regions at 90%, 95%, 99%, and 3cr C.L. for 
2 d.o.f. The left panel corresponds to Fig. [6j whereas in the right panel the new 
background from 13 C has been taken into account. 



finds 

Aw?, 2 L 

sin 2 29 12 sin 2 — ^- ~ 1 - P ee - 2P ee sin 2 9 13 (A.2) 

For a given survival probability Eq. (1A.2|) implies a negative correlation between sin 2 13 
and sin 2 2^12 . This trend is visible also in the actual fit, see Fig. \A2\ where the allowed 
region for these two parameters is shown for fixed Am 2 . One observes that increasing 
sin 2 9 13 can be compensated to some extent by decreasing sin 2 12 (within the "light 
side" sin 2 6*12 < 0.5). This, however, is in disagreement with solar data, which provide a 
stable lower bound on sin 2 9 12 due to the fundamentally different conversion mechanism 
(MSW matter effect in the sun versus vacuum oscillations). The combined analysis of 
KamLAND and solar data leads essentially to the intersection of the two individual 
allowed regions, which further improves the bound on sin 2 13 . 



Appendix B. Implications of a new background in KamLAND 

After this paper has been published it was realised that a background from the reaction 
13 C(a, n) 16 contributes to the KamLAND data, which has not been taken into account 
in the first version of Ref . 15) , on which our analysis is based. In this appendix we show 
that the impact of this new background on our results is small, and we present the 
updated results, where the changes are noticeable. 



Status of global fits to neutrino oscillations 37 




Figure B2. Left panel: A\ 2 from KamLAND data allowing for negative values of 
sin 2 6*i3 with and without accounting for the 13 C background. Right panel: A\ 2 as 
a function of sin 2 #13 for various data samples. The 13 C background is (is not) taken 
into account for the solid (dashed) curves. 



The reaction 13 C(a, n) 16 leads to 10.3±7.1 events above the 2.6 MeV threshold in 
KamLAND, and hence the total background is increased to 17.8 ± 7.3 events. This new 
background is mainly concentrated around 6 MeV, and the main effect of subtracting 
these events from the reactor data is that the relatively high data point at 6 MeV 
visible in Fig. IA1I is moved from a value of 1.1 to 0.8, in better agreement with the 
prediction for oscillation. Indeed, we now find xiun = 9-5 instead of 13.5, i.e. the 
quality of the fit improves. The impact of the new background on the determination 
of the solar parameters is shown in Fig. IB1I The best fit point for Am 2 0L moves from 
8.1 x 10~ 5 eV 2 to 7.9 x 10~ 5 eV 2 . However, the allowed region around the best fit 
point is very stable. The fact that the best fit point occurs now for sin 2 # S0L > 0.5 
has no statistical significance. The low-Am 2 OL solution is slightly more disfavoured, 
with Ax 2 = 8.6 instead of 7.4, whereas the the high-Am 2 OL solution becomes somewhat 
better, with A% 2 = 8.8 instead of 11.2. Since these solutions are ruled out anyway by 
solar data, the overall changes in the global fit are very small. 

The impact of the new background on the bound on #13 is illustrated in Fig. IB2I The 
left panel shows the \ 2 °f KamLAND data without imposing the constraint sin 2 6* 13 > 0. 
We observe that the best fit occurs now closer to the physical region, and the A% 2 
between sin 2 # 13 = and the best fit point decreases from 3.5 to 1.75. This again 
shows that the quality of the fit improves. As a consequence the bound on sin 2 #13 from 
KamLAND (constraining the fit to the physical region) becomes slightly weaker. From 
the right panel one can see that the impact on the bound from global data is rather 
small. Finally, in Tab. IB II we update the results of the global three-neutrino analysis 
given in Tab. [IJ taking into account the new background in KamLAND data. Only the 
numbers for Am 2 0L and sin 2 #13 change slightly. 
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parameter 


best fit 


2a 


3a 


4a 


Am§! [10" 5 eV 2 ] 


7.9 


7.3-8.5 


7.1-8.9 


6.8-9.3 


Amli [10- 3 eV 2 ] 


2.2 


1.7-2.9 


1.4-3.3 


1.1-3.7 


sin 2 


0.30 


0.25-0.34 


0.23-0.38 


0.21-0.41 


sin 2 #23 


0.50 


0.38-0.64 


0.34-0.68 


0.30-0.72 


sin 2 0i3 


0.000 


< 0.031 


< 0.051 


< 0.073 



Table Bl. Updated version of Tab. Q]taking into account the background from 13 C 
in KamLAND: Best-fit values, 2cr, 3a, and 4u intervals (1 d.o.f.) for the three-flavour 
neutrino oscillation parameters from global data including solar, atmospheric, reactor 
(KamLAND and CHOOZ) and accelerator (K2K) experiments. 



Appendix C. June— 2006 update 



We update our three-neutrino oscillation parameter analysis using the most recent 
data: the first MINOS result s us ing the NuMI Beam w hich have been presented at 



the Neutrino2006 conference 



141 



collaboration 145 



the last K2K results 



142 



the new update of the 



Standard Sol ar M odel (SSM) |143l . 1 1441 ] and also latest measurements from the SNO 



First MINOS data 

MINOS is a long-baseline experiment that searches for disappearance in a 
neutrino beam with a mean energy of 3 GeV produced at Fermilab. It consists of a near 
detector, located at 1 km from the neutrino source and a far detector located at the 
Soudan Mine, at 73 5 km from Fermilab. Recently data corresponding to 1.27x 10 20 p.o.t. 
have been released 



141 



slightly more than the final K2K data sample. In the absence 
of oscillations 239 ± 17 events with E < 10 GeV are expected, whereas 122 have 
been observed, which provides a 5.9a evidence for disappearance. In our re-analysis we 
use spectral data divided into 15 bins in reconstructed neutrino energy, and o ur a llowed 
region from MINOS-only is in very good agreement with the official result 141] . The 



values of the oscillation parameters from MINOS are consistent with the ones from K2K, 
as well as from SK atmospheric data. 

Last K2K data 

Recently the K2K collaborati on h as published details of the analysis of their full 



data sample (K2K-I and K2K-II) |l42j |. The data have been taken in the period from 
June 1999 to November 2004 and correspond to 0.922 x 10 20 p.o.t. Without oscillations 
15811 J events 

are expected whereas only 112 events have been observed. Out of these, 58 
events are single-ring events where the reconstruction of the neutrino energy is possible, 
we use these events to perform a spectral analysis of the K2K data, as described in 
Sec.E21 



New SNO-salt data and SSM update 
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We have also updated our solar data analysis in view of the recently released update 
of the SSM and the rece nt SNO d ata in its salt phase. Taking into account new radiative 



opacities, Bahcall et al. |143l . Il44| obtain new solar neutrino fluxes, neutrino production 
distributions and solar density profile. Among the different solar models they present, we 
have adopted the one denoted as BS05(OP), since it agrees better with the helioseismic 
measurements. For completeness we have checked also the alternative model choice 
BS05(AGS,OP) obtained from recently determined heavy element abundances. The 
change to this model, however, has a negligible effect on the result for the oscillation 
parameters found for the preferred BS05(OP) option. (For our analysis the main 
difference between these two models is the predicted boron flux, which is now more 
precisely determined by SNO data). 

In addition to the new SSM we include in the pr esen t analysis the latest SNO 



neutrino flux (CC, NC, ES) measurements given in Ref. [145] based on 391 days of data 
for the SNO phase II. With respect to the previous result (254 days) the value of the 
CC to NC flux ratio changed from <Pcc/<Pnc = 0.306 ± 0.026 ± 0.024 to the slightly 
higher value <Pcc/<Pnc = 0.340 ± 0.023 ± 0.030. 

Updated three— neutrino analysis 

The results of our updated global three-neutrino analysis are summarized in Fig. IC1I 
and Tab. IC1I Here we simply highlight the main differences with respect to our previous 
results: 

• An increase in the best fit value for the "atmospheric" mass splitting Am^ and a 
narrower allowed range for this parameter due to the inclusion of the first MINOS 
data which mainly lead to a tighter lower bound (see Fig. IClj) . The allowed range 
of #23 is dominated by atmospheric data as before, since MINOS provides a very 
weak constraint on the mixing angle, similar to K2K. 

• A small shift in the allowed interval for the "solar mixing angle" sin 2 # 12 , due to the 
larger value of last SNO-salt fluxes compared to its previous measurements and the 
change in the SSM. The so-called high-LMA solution around Am 21 ~2x 10~ 4 eV 2 
is now slightly less disfavored, with A% 2 « 15. This follows from the somewhat 
higher value of (ficc/^NC, as becomes clear from the contours of (pcc/^NC shown 
in Fig. [5j 

• A slight improvement in the #13 bound, mainly because of the larger value of Arn^j 
preferred by MINOS data. In the global analysis also the updates in the solar fit 
contribute to the improved bound. 



The small parameters: a and #13 

In Fig. IC2l we show the \ 2 as a function of the a and a sin 26u parameters relevant 
for three flavour effects in future long-baseline experiments, where a denotes the ratio 
of "solar" to "atmospheric" mass-squared differences. With the new data we obtain the 
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Figure CI. Projections of the allowed regions from the global oscillation data at 90%, 
95%, 99%, and 3<r C.L. for 2 d.o.f. for various parameter combinations. Also shown is 
Ax 2 as a function of the oscillation parameters sin 2 #12, sin 2 #23, sin 2 #13, Am 2 !, Am 2 ^, 
minimized with respect to all undisplayed parameters. Dashed lines and empty regions 
correspond to the global analysys before this update, while solid lines and coloured 
regions show our most recent results. 



parameter 


best fit 


2a 


3a 


4a 


Am 2 21 [10" 5 eV 2 ] 


7.9 


7.3-8.5 


7.1-8.9 


6.8-9.3 


Amli [l(T 3 eV 2 ] 


2.6 


2.2-3.0 


2.0-3.2 


1.8-3.5 


sin 2 9i2 


0.30 


0.26-0.36 


0.24-0.40 


0.22-0.44 


sin 2 623 


0.50 


0.38-0.63 


0.34-0.68 


0.31-0.71 


sin 2 #13 


0.000 


< 0.025 


< 0.040 


< 0.058 



Table CI. 2006 updated version of Table 1. Best-fit values, 2er, 3cr, and 4cr 
intervals (1 d.o.f.) for the three-flavour neutrino oscillation parameters from global 
data including solar, atmospheric, reactor (KamLAND and CHOOZ) and accelerator 
(K2K and MINOS) experiments. 



following best fit values and 3a allowed ranges: 

a = 0.030, 0.024 < a < 0.040, (C.l) 

a sin 2#i2 = 0.028 , 0.022 < a sin 26 12 < 0.037 . 

The local minimum at A% 2 w 15 visible in Fig. IC2I appears because of the somewhat 
weaker rejection of the high-LMA solution due to the new SNO CC/NC value. 
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Figure C2. Left: A% 2 from global oscillation data as a function of a = Am 2 OL /Am^ TM 
and a sin 2^i2 after the inclusion of MINOS data. Right: 90% C.L. upper bound on 
sin 2 (2 d.o.f.) from the combination of all neutrino oscillation data as a function of 



Am 2 TM . 



The right panel of Fig. IC2I illustrates how the bound on 9 ri emerges from an 
interplay of the global data. In particular, the role of solar and KamLAND data is 
clearly visible, even after the improvement of the lower bound on Am^ from MINOS. 
We find the following bounds at 90% C.L. (3a) for 1 d.o.f.: 

0.033 (0.071) (solar+KamLAND) 
sin 2 #i 3 <<{ 0.026 (0.054) (CHOOZ+atmospheric+K2K+MINOS) (C.2) 
0.020 (0.040) (global data) 

Appendix D. September— 2007 update 



We have update d ou r analysis including the new data released by the MINOS |l46 



and KamLAND [147] collaborations. New MINOS data have been collected from June 
2006 to July 2007 (Run-IIa), and they have been analyzed together with the first data 
sample (Run-I), with a total exposure of 2.5xl0 20 p.o.t. In total, 563 events have 
been observed at the far detector, while 738±30 events were expected for no oscillation. 



; exp 
1471 



The most recent data from the KamLAND experiment [1471 ] correspond to a total 
exposure of 2881 ton-year, almost 4 times larger than 2004 data. They provide a 
very precise measurement of the solar neutrino oscillation parameters, mainly the mass 
squared splitting. Apart from the increased statistics this is also due to the reduction 
of systematic uncertainties. Thanks to the full volume calibration the error on the 
fiducial mass has been reduced from 4.7% to 1.8%. The main limitation for the Am^ 
measurement comes now from the uncertainty on the energy scale of 1.5%. The analysis 
of the new data is performed in a similar way as described in appendix A. We use 
the data binned in equal bins in 1/E to make optimal use of spectral information. As 
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Figure Dl. Left: Allowed region in the (sin 2 623, Am^) plane before (lines) and after 
(coloured regions) the inclusion of the new MINOS data. Right: Allowed region in the 
(sin 2 0i2, Ato 2 ^) plane before (lines) and after (coloured regions) the inclusion of the 
new KamLAND data. 



previously we restrict the analysis to the prompt energy range above 2.6 MeV to avoid 
large contributions from geo-neutrinos and backgrounds. In that energy range 1549 
reactor neutrino events and a background of 63 events are expected without oscillations, 
whereas the observed number of events is 985. 



The Borexino collaboration has also presented their first data 148]. While they 



provide the first real time detection of Berilium-7 solar neutrinos, and an important 
confirmation of the standard solar model and the large mixing oscillations, they currently 
do not affect the determination of neutrino oscillation parameters. 

Our results for the atmospheric and solar neutrino oscillation parameters are 
summarized in Fig iDll In the left panel we show the allowed region and A% 2 profiles 
for the atmospheric parameters (sin 2 #23, Am^) before and after the inclusion of 
the MINOS-2007 data. One can appreciate a decrease in the best fit value and 
a narrower allowed range for the "atmospheric" mass square splitting Am^. This 
illustrates the increasing role of long baseline accelerator data in the determination of 
the "atmospheric" splitting, a trend that will get more pronounced in the future. The 
right panel gives the analogous plot for the solar neutrino oscillation parameters (sin 2 # 12 , 
AmU, where one can see how the 2007 KamLAND data provides an improvement in 
the determination of the "solar" mass square splitting as well as the lower bound on 
sin 2 6 l 1 2, while the upper bound is still dominated by SNO solar neutrino data. 

The information concerning the mixing angle #13 is shown in Fig JD2l In the left 
panel we have plotted the A% 2 profile as a function of sin 2 # 13 from the analysis of solar 
+ KamLAND, atmospheric + K2K + MINOS + CHOOZ, and also from the global 
analysis of all the data samples. In the global analysis we find a slight weakening of 
the upper bound on sin 2 #13 from 0.04 to 0.05 at 3a. The reason for this is two-fold. 
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Figure D2. Left: A% 2 profiles as a function of #13 from the analysis of different data 
samples. Right: 90% C.L. upper bound on sin 2 # 13 (2 d.o.f.) from the combination of 
all neutrino oscillation data as a function of Am 2 .,. 



parameter 


best fit 


2a 


3a 


Am 2 21 [10" 5 eV 2 ] 


7.6 


7.3-8.1 


7.1-8.3 


Arr4 [l(T 3 eV 2 ] 


2.4 


2.1-2.7 


2.0-2.8 


sin 2 9i2 


0.32 


0.28-0.37 


0.26-0.40 


sin 2 6> 23 


0.50 


0.38-0.63 


0.34-0.67 


sin 2 6»i3 


0.007 


< 0.033 


< 0.050 



Table Dl. 2007 updated version of Table 1. Best-fit values, 2a and 3<r intervals 
(1 d.o.f.) for the three-flavour neutrino oscillation parameters from global data 
including solar, atmospheric, reactor (KamLAND and CHOOZ) and accelerator (K2K 
and MINOS) experiments. 



First, the shift of the allowed range for Am 2 x to lower values implies a slightly weaker 
constraint on sin 2 #13, and second, the combination of solar and KamLAND data prefers 
a slighlty non-zero value of sin 2 9 13 which, though not statistically significant, also results 
in a weaker constraint in the global fit. In the right panel we show the 90% C.L. upper 
bound on sin 2 #13 from the combination of all data samples. No significant improvement 
has been otained here, besides the narrower range for Am^. 

A summary of the updated neutrino oscillation parameters is given in Tab. ID It 
Including the 2007 data we find the following best fit values and 3a allowed ranges for 
the parameters characterizing three flavour effects in future long-baseline experiments, 
namely #13, a and a sin 29i 2 : 

{ 0.051 (0.084) (solar+KamLAND) 
0.028 (0.059) (CHOOZ+atmospheric+K2K+MINOS) (D.l) 
0.028 (0.050) (global data) 

a = 0.032, 0.027 <a< 0.040, (D.2) 

a sin 29 12 = 0.029 , 0.024 < a sin 29 12 < 0.037 . 
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